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Abstract 

Berezinskii-Kosterlitz-Thouless (BKT) transition is one of the instability mech- 
anisms for the Tomonaga-Luttinger liquid. But in the BKT transition, there are 
logarithmic-correction problems, which make it very difficult to treat BKT transi- 
tions numerically. We propose a method, "level spectroscopy", to overcome such 
difficulties, based on the renormalization group analysis and the symmetry consid- 
eration. 



1 Introduction 

Tomonaga-Luttinger (TL) liquid [p], @] is an important concept for one dimensional (ID) 
quantum systems (spin, electron systems, nanotube, etc. at T = 0) with the continuous 
symmetry (U(1),SU(2), etc.), and TL liquid is closely related to several 2D classical models 
(6-vertex model, classical XY spin, superconducting film, roughening transition, etc. at 
T^O). Therefore, it becomes crucial to understand instabilities of TL liquid. 

One of the instabilities of TL liquid is the Berezinskii-Kosterlitz-Thouless (BKT) tran- 
sition [[| |], m . Although the BKT transition has been well known, it has not been recog- 
nized that there are pathological problems to analyze numerically BKT transitions. One 
of these problems is that the finite size scaling technique || , which is successful for second 
order transitions, can not be applied for BKT transitions |7], [|, since there are logarith- 
mic corrections from the marginal coupling. Recently, combining the renormalization 
group calculation and the symmetry, we have developed a remedy, "level spectroscopy" , 
to overcome these difficulties. 

In the next section, we compare the BKT transition with the second order transition. 
In §3, we introduce the concept of the level spectroscopy and how to use it. In §4, we 
deduce the level spectroscopy from the sine-Gordon model, which is an effective model 
to describe the BKT transition. For readers who are not familiar with the sine-Gordon 
model, we recommend to read §5 (physical examples) before §4. 
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2 BKT transition versus second order transition 

First, we briefly review the renormalization concept. Let us consider a d- dimensional 
classical Hamiltonian (or an action for a (d — 1) + 1 dimensional quantum system) |9| 

3 

where Ho is scale invariant, ipj(x) is a local order parameter, yj is an effective coupling 
constant (or some external field), a is a short-range cutoff. 

When changing a scale as a' = ba, effective local order parameters change as ipj = 
b x iipj (xj\ scaling dimension) and according to this, effective couplings change as y'j = 
b d ~ Xj yj. If yj = 0, then y'- remain (fixed point). We call the case y'- diverging for b — > oo 
(i.e. Xj < d) as relevant, whereas the case y'j converging (i.e. Xj > a) as irrelevant. 

At the fixed point, the system is scale invariant, therefore the correlation length is 
infinite £ = oo. Although regions which flow into the fixed point are not strictly scale 
invariant, £ = oo in these regions (critical regions). 

For the infinitesimal scale transformation a' = ae dl w a(l + dl), one can treat eq. (P 
perturbatively. 



2.1 Second order transition 

The scaling equations for second order transitions are 

< ^p- = (d- Xj ) yj (l) =(3^(1)}) {l = \n{L),L: system size). (2) 

As an example, let us consider the two scaling field case. When ipi is a relevant operator 
and tp2 is an irrelevant one, the renormalization group flow is given in Fig. [T] a. 

Since the second order transition occurs where the sign of the j3 function (for a relevant 
coupling 2/i ) changes, one can use this fact to determine the critical point 0. 




Figure 1: Renormalization flows for (a) second order transition, (b) BKT transition. 
Shaded parts are critical region. 
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2.2 BKT transition 



The BKT transition (d = 2) is described by the following RG equation |J 



dyijl) 
dl 



dl 



Vl{l)V2(l) , 



(3) 



where fixed points are y 2 = 0. All of the points in \y 2 \ < y\ are renormalized to y 2 = 0, 
thus in this region correlation lengths are infinite (£ = oo) or massless. Other regions are 
massive (£ finite), except y 2 = (see Fig [I] b). We call the region within \y 2 \ < y\ as the 
massless region, the lines y 2 = ±j/i as the BKT critical lines, y\ = as the Gaussian fixed 
line. The point 2/1 = 2/2 = has a special meaning (two BKT lines and one Gaussian line 
intersect), and we call it as the BKT multicritical point. 

Note that on the BKT line, the coupling y\ is marginal (i.e. X\ = 2) and it behaves as 
2/i (Z) = 1/7 = 1/ ln(L). Another important fact is that on the Gaussian fixed line (y 2 = 0), 
the scaling dimension of ip 2 is varying with y\. 

2.3 Comparison between second order and BKT transition 

1. Critical region (£ = 00) is isolated in the second order transition, whereas it is 
extended in the BKT transition. 

2. Zero point of the f3 function corresponds to the critical point in the second order 
transition, whereas for the BKT, it has no special meaning. 

3. Since the renormalization group behavior on BKT critical lines (y 2 = ±2/1) is 
marginal, there appear logarithmic corrections 1/ ln(L), e.g., in the correlation func- 



therefore finite size effects are very large. 

All of them make it very difficult to treat the BKT transition numerically (N.B. 
universal jump is also affected by logarithmic corrections, but 0(l/(lnL) 2 )). 

3 Level spectroscopy 

The BKT transition occurs where some physical quantities change from irrelevant to 
relevant. Thus it is useful to investigate scaling dimensions near marginal. In fact, on 
the Gaussian line (i.e., without interaction), several scaling dimensions cross at the BKT 
multicritical point (see Fig. |2] a, or eqs. (|TTD, ([TED). 




(4) 



or in the energy gap for the finite size system as |10|, [TT 
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Based on the conformal field theory (CFT) in 2D [I2"|, one can obtain the scaling 
dimensions using the energy gap for the finite system. The relation between the energy 
gap AEj of the system size L and the scaling dimension Xj is 13] 



2TTVXj 



(5) 



where v corresponds with a renormalized Fermi velocity, or a spin wave velocity. 

Each excitation can be classified by quantum numbers (m (magnetization or electron 
density, related with U(l) symmetry), P (parity), q (wave number)). 

In the normal BKT transition, there is no symmetry breaking in the massive phase. 
But in general, the BKT transition may be combined with a discrete symmetry, and there 
occurs a Z p discrete symmetry breaking in the massive phase. 

Procedure to use level spectroscopy 

1 . Classification of BKT transitions with the discrete symmetry 

(a) BKT transition without symmetry breaking (§3.1, §3.3) 

(e.g. 2D classical XY, integer S XXZ quantum spin chain ) 

From table [l], choose the excitation with some quantum numbers, then from 
the level crossing, determine BKT transition line. 

(b) BKT transition with the Z 2 symmetry breaking (§3.2) 

(e.g. 6 vertex model, half-integer S XXZ quantum spin chain ) 

From table |2], choose the excitation with some quantum numbers, then from 

the level crossing, determine BKT transition line. 

(c) BKT transition with the Z p (p > 2 integer) symmetry breaking 
(e.g. 2D p-state clock model) 

Although we do not explain the Z p case here, this case has been discussed in 



14 1 . Note that there is a difference for p even or odd case. 



2. Checking of the universality class (elimination of logarithmic corrections) 

(a) Scaling dimension (§3.4) 

From tables [I], @> choose the excitations eliminating logarithmic corrections 
each other, and check the universality class. 

(b) Central charge c = 1 (§3.5) 



3.1 Normal BKT transition 

In the normal BKT transition, there is no symmetry breaking. In table [I], we show the 
relation between the scaling dimension (or excitation (|5|)) and quantum numbers fUl . In 
the neighborhood of the scaling dimension x = 2, on the BKT line 2/2(0 = ±Vi(l), there 
is a level crossing of excitations with quantum numbers (m = ±4, P = l,q = 0) and 



French-Japanese Symposium in Fukuoka, 2001 



5 



(m = 0, P = 1, q = 0) (see table [I], Fig. |2|), thus we can determine the BKT transition 
line. 

Next, since the ratios of logarithmic corrections (1/lnL) on the BKT transition line 
are 2 : 1 : — 1 (corresponding to xq )COS) xo.sm, ^±4,0 m table 0), we can eliminate logarithmic 
corrections and check the universality class. 

This level crossing on the BKT transition reflects the (hidden) SU(2) symmetry p0| , 
thus it is correct up to higher order loops. We can see this SU(2) symmetry explicitly, 
including the twisted boundary condition (see §3.3). 




Figure 2: Scaling dimension x in the neighborhood of BKT transition (a) y 2 = (on the 
Gaussian fixed line) (b) y 2 ^ 0. 



m P q 


BC 


X 


operator in s.G. 


abbr. 


±2 1 
-1* 
1* 


PBC 
TBC 
TBC 


1/2 -y x (/)/4 

1/2 + ^(0/4-^(0/2 

1/2 + ^(0/4 + 1/2(0/2 


sin(0/V2) 
cos(0/ v / 2) 


^±2,0 


±4 1 
1 
0-10 
1 


PBC 
PBC 
PBC 
PBC 


2-2/i(0 

2- yi (0(l + 4t/3) 
2 + Vl {l) 

2 + 2j/i(0(l + 2*/3) 


exp(±i4 v / 2^) 
marginal 
sin (720) 
cos(\/20) 


^±4,0 
•Emarg 

Xq,cos 



Table 1: Renormalized scaling dimensions x and quantum numbers (m, P, q) for the nor- 
mal BKT transition. PBC denotes periodic boundary condition, TBC twisted boundary 
condition, t is a distance from the BKT critical line, defined as 2/2 (0 — ±2/i(0(l + t)- 
Couplings yi, i/2 follow renormalization group equations @, and on the BKT transition 
\me yi (l) = ±y 2 (l) = l/\n(L/L ). 

3.2 BKT transition with Z2 symmetry breaking 

Next we consider the BKT transition coupled with a discrete symmetry. Especially in the 
BKT transition with the Z 2 symmetry breaking (in the massive region), the level crossing 



French-Japanese Symposium in Fukuoka, 2001 



6 



of the lowest excitations in each region corresponds to the phase boundary. 

We introduce a quantum number corresponding to the Z2 symmetry (wave number 
q — 0, 7r etc.). In table ^| we summarize quantum numbers and excitations [15, 16 . 
Level crossings on BKT transition lines can be observed not only x = 2 (N. B. quantum 
numbers differ from table [I]), but also x = 1/2, where the excitations with quantum 
numbers (m = 0, q = 7T, P = ±1) and (m = ±1) show a level crossing (Fig. |3| a, table |2|). 

About the Gaussian fixed line 2/2 = 0, each of the Z2 symmetry broken massive phases 
has different parity, thus it occurs a level crossing of excitations m = 0, q = n, P = ±1 
(Fig. |b). 



m P q 


X 


operator in s.G. 


abbr. 


±1 

-1 i 

1 7T 


l/2- yi (l)/A 

l/2 + yi {l)/i-y 2 (l)/2 

l/2 + yi(0/4 + y2(0/2 


exp(±iV29) 

sin( v / 20) 

cos(\/20) 


£±1,0 

•E0,cos 


±2 1 
1 
0-10 
1 


2-2/i(0 

2-^(0(1 + 44/3) 
2+2/i(0 

2 + 2^(0(1 + 24/3) 


exp(±i2V2#) 
marginal 
sin(2V2» 
cos(2 v / 20) 


X±2,0 
•Emarg 

%0,cos2 



Table 2: Renormalized scaling dimensions x and quantum numbers (m, P, q) for BKT 
transition with Z 2 symmetry. 




3.3 Twisted boundary condition (TBC) 

In the previous section case, one can distinguish two massive phases on the both side of 
the Gaussian fixed line, by using Z 2 symmetry. How to distinguish two massive phases 
in the normal BKT transition? Using the twisted boundary condition (TBC) method 
171 115), we can clarify the hidden Z 2 x Z 2 symmetry. 
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TBC is expressed by the sine-Gordon language as (see §4) 

exp(±iV26(z , z\ + L)) = — exp{±iV26(z , Zl )), (6) 

or in the quantum spin language (S=l case, see §5) 

q x ,y — q x >y qz — qz (>j\ 
°L+j — °j > °L+j — °j ■ \ ' ) 

We can also define the discrete (inversion) symmetry P* : Sj <-> Sl+i-j under TBC 
(N. B. differs from the parity under PBC). Corresponding to this quantum number, we 
can observe the level crossing of m = states under TBC at the Gaussian fixed line y% = 
(see table |l|). 

Besides, using the TBC method, we can clarify the hidden SU(2) symmetry on the 



BKT critical line [|T9| , |20| . That is, we can determine BKT line by the level crossing 
between the excitations under PBC q = 0, m = ±2 and that under TBC m — 0. (see 
table 0). 

3.4 Universality class (scaling dimension) 

Finally, in order to check the consistency of our method, we should eliminate logarith- 
mic corrections from scaling dimensions (critical indexes). There are several methods to 
eliminate logarithmic size corrections on BKT lines. 
For the BKT transition with Z2, relations 

(x ,sin + X0,cos) x £±1,0 = 1/2, (8) 
x±2,o/x±ifl = 4, (9) 

are correct up to 0(l/(lnL) 2 ), and they also apply all over the critical region. 

Similarly, for the normal BKT case, combining TBC, we can check the universality 



class as the above method 18 . 



3.5 Central charge 

The BKT critical region (massless region £ = 00) also can be characterized by the central 
charge c = 1. Numerically the central charge c is obtained from the ground state energy 
for the finite system as J2"I| 

E g {L) = e 9 L - ™ (10) 

(N. B. the central charge is also obtained experimentally from the specific heat [p^] .) 

Although there are logarithmic corrections in the effective central charge c obtained 
from eq. (|10|), they are small enough 0(l/(lnL) 3 ) 0], thus we can neglect them. 

Note that the effective central charge, obtained as eq. flip]) , changes rapidly from c = 1 



(massless) to c = (massive) 23, 24 . 
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4 Sine-Gordon model 

There are several effective models to describe the BKT transition. 

f . One of the effective models for BKT transitions is a sine-Gordon model, 

Z = Jv4>exp(- Jd 2 xC), £ = ^(V0) 2 + ^cospv^, (11) 

where p is integer, K is related with y x = 2(Kp 2 /4 — 1) in eq. (|3|). In addition to 
eq. (p!l|), we require following features: 

(a) Compactification: = + \^2n. 

(b) Canonical field 9 to 

<9 X .0 = -d y (iK0), d y <P = d x (iK0), (12) 

(c) Compactification of 9: 9 = 9 + \f2~-K. 

This represents the continuous U(l) symmetry, because eq. (|11]) is invariant 
under 9 — > 9 + const. 

(d) Discrete symmetry 

Eq. (|TT1) is invariant under the change 0^0 + V2ir/p, which corresponds to 
the discrete Z p symmetry. 

2. Free case (y 2 = 0) @ 

(a) Correlation functions for and 9 

It is convenient to describe coordinates in 2D as complex variables: z = x + iy. 
Then, correlation functions for 0, 9 are 

mm) = -| ln (M), w , W0)) = __L hl (M), 

M*)0(O)> = -^argz. (13) 

(b) Vertex operators 

O m<n = exp(i^/2m9) exp(iv2n(j)) (m, n : integer). (14) 

(c) Marginal operators 

Omarg = (9,0) 2 + (9,0) 2 . (15) 

(d) Correlations and scaling dimensions for vertex operators 



(O min (z)0_ m _„(0)) = exp 



-2x m , n In f ^- j - 2z/ mi „(arg z + tt/2) 



(16) 
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Therefore, scaling dimensions are given by 

(17) 



1 (m? \ 



thus x m)n are varying with coupling K, whereas Z m>n (relating with the wave 
number in the ID quantum system) are fixed. 

(e) Correlation and scaling dimension for marginal operator 

(O m arg(z)O marg (0)) OC \z\~ 4 . (18) 



Thus, the scaling dimension for the marginal operator is x marg = 2, I = 
3. Interacting case (y 2 7^ 0) 

(a) Parity {(f) — ► —0; corresponding to the space inversion in ID quantum system) 
Considering parity, we choose operators cos(\/2n0) and sm(v2?T<0) ■ 

(b) Renormalization from t/j term 

Couplings yi and y 2 are renormalized as eqs. (§). This affects all the scaling 
dimensions, according as eq. fllT|). 



4.1 Normal BKT (p = 1) 

In this case, the yi coupling term in (|ll|) becomes relevant at K = 4 on y% = 0. 

1. Level crossing at K = 4, y 2 = 

On Gaussian line at K — 4, according to eqs. ([H]), (|18l), five operators ((m = 
±4,P = l,g = 0),(m = 0,P = 1, g = 0), (m = 0, P = -1, q = 0), (m = 0, P = 
1, g = 0)) have the same scaling dimensions (x = 2). 

2. Hybridization (y 2 7^ 0) 

The operator cos V2(J) and the marginal operator are affected from the renormal- 
ization of 2/2, since they have the same symmetry (m, P), they hybridize each other 
by the 2/2 term |14j (calculation can be more simplified with the operator product 
expansion (OPE) [pl|). Combining the renormalization from yi, there remains a 
level crossing on the BKT lines (see table |]). 



4.2 BKT with Z 2 (p = 2) 

In this case, the y 2 coupling term in (|11|) becomes relevant at K = 1 on y 2 = 0. 

1 . Level crossing at K = 1 , y 2 = 

On Gaussian line at K = 1, besides the level crossing at x = 2, the four operators 
((m = ±1,P = —l,q = ir),(m = 0,P = ±l,g = 7r)) have the same scaling 
dimensions (x = 1/2), 
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2. Level split (y 2 ^ 0) 

With the cos 2\^2(j) coupling, scaling dimensions for cos \/2(j) and sin \/2(f) are s P n ^ 



15, [L6fl. Combining the renormalization from yi, there remains a level crossing on 



the BKT lines (see table 



5 Physical examples 

Here we show physical examples for the quantum spin and the electron chains. Note 
that spin systems obey the commutation relations, whereas electron systems obey the 
anti-commutation relations. Thus, for the spin chain, one can directly relate the quantum 
numbers to those of the sine-Gordon model, whereas for the electron case, we should 
choose an appropriate boundary condition according to evenness or oddness of quantum 
numbers (selection rule). 

5.1 S=l spin chain 

First we consider the S=l bond-alternating XXZ chain, 



L 



# = 8(-lY)(SfS? +1 + S?S? +1 + AS*S* +1 ). (19) 

i=i 

In this case, quantum numbers are defined as the magnetization m = Sj , the parity P 
for the space inversion Sj <-> Si+i-j, the wave number q for the translation by two sites 
Sj — > Sj + 2- The level crossing of excitations (m = ±4, P — l,q — 0) and (m = 0, P — 
l,q = 0) corresponds to the BKT phase boundary. The detailed analysis of the phase 
diagram and the universality class for this model is given in |26], Note that using 



TBC method, one can improve the accuracy [17, 19 



5.2 S=l/2 spin chain 

Next we consider the S=l/2 XXZ chain with the next-nearest-neighbor interaction. 

H = X>i,i+i + "/'„.•>). K,i = S*S* + SfS] + AS*S*. (20) 

3=0 

In this case, the level crossing of excitations (m = ±1, q = tt) and (m = 0, P — ±1, q — it) 
corresponds to the BKT phase boundary. The detailed analysis of the phase diagram and 



the universality class for this model is given in [IB . 



5.3 Electron system: selection rule 

We briefly review the history of the selection rule and the boundary condition of the ID 



fermion model. Using the Jordan- Wigner (non-local) transformation, Lieb et al. [27 
have studied the exact mapping from a S=l/2 spin chain, which is equivalent to a hard 
core boson, to a spinless fermion chain. They have pointed out that according to the 
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oddness or evenness of the fermion number, one should use the PBC or TBC In the sine- 
Gordon model (phase Hamiltonian) mapped from the spinless fermion, or the bosonization 
language, Haldane has written a systematic review. In that paper, he has introduced 
a new quantum number, current, J = Nl — N R , and he has written a selection rule for 
fermion numbers and current number. One can see from eq. (3.54) in ||, that boundary 
condition should change according to these quantum numbers (although in only the 
forward scattering case was discussed, it is possible to include the umklapp interacting 
case). The extension from the spinless fermion case to the electron chain (fermion with 
the spin freedom) is straightforward, considering two species of spinless fermion chains. 
From the another point of view, the selection rule has been found in the field of Bethe 
Ansatz. For the Hubbard model, it is described by Woynarovich . 



Returning to the concrete procedure, we consider the spinless fermion case given by 



L L 



H = -t X, [ c i c i+i + c j+i c i\ +V z2 r ] r j r j-\<\, :- ( 21 ) 

3=1 3=1 

For the ground state, we choose the following boundary condition. When the particle 
number is odd, we assume PBC. But when N is even, the ground state is two-fold de- 
generate. In order to remove this degeneracy, we assume TBC Cj +L = —Cj. We write 
the fermion number for the left mover as Nl, and for the right mover as N R (we do not 
include q = fermion as Nl or Nr). We define another numbers m and n as 

Nr - Nr 

m = N L + N R -N , w= L 2 R . (22) 

When the particle number is odd, N is the particle number of the ground state minus 1 
with PBC, and when the particle number is even, N is the particle number of the ground 
state with TBC. m means the change of the particle number. These numbers relate to 
the boundary condition of the phase field [0] 

<f>(x + L) = <j)(x) - V2nm, 9(x + L) = 6(x) - V2nn. (23) 

The low-lying excitation spectrum of the system is given by 

E mtn (h L , h R ) - E = —x min + ~j-(n L + n R ) (24) 

where E is the ground state energy, v is the sound velocity, and x mtn is given by eq. ( |T7D 
which is the scaling dimension of the operator (0). The second term of eq. (]24]) gives the 
sound wave collective excitation and ul,r are non-negative integers. The wave number of 
this excitation is given by 

/ 2lT \ 27T 

q = - ( 2k F + — mj n - — (n L - n R ), (25) 

where hp is the Fermi wave number. Since Nl and N R are integer, the number n ± m/2 
must be so. Thus when m is an even integer, n is an integer, and when m is an odd 
integer, n is half odd integer ((— l) m = (— l) 2n ). For Tomonaga and Luttinger models, 
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Nl and N R are good quantum numbers. When some interactions which do not conserve 
Nl and N R , such as the umklapp scattering, are introduced, m remains a good quantum 
number but 2n does not. Although only the parity (— l) 2 ™ is conserved, it does not violate 
the selection rule mentioned above. The instability of the BKT transition for the model 
fl2ip occurs at V = It for the half filling case J2jC^Cj = L/2, and this instability stems 
from the umklapp scattering. 

Next we consider the electron case with spin freedom, such as the ID extended Hub- 
bard model 

L L 

H = ~tJ2Yl [ c i<T C *+l» + c i+l«rC*r] + U n j] n jl 
3=1 <x=U ' 3=1 

L 

+ v Yl Kt + n 3i) K+iT + n 3+n) » ( 26 ) 

3=1 

where nj a = Cj a Cj a . In this case, we have four particle numbers of the left and the right 
movers with spin er =fj, Nia and Nru . We assume PBC when the particle number of 
the ground state is AM + 2 (M is integer), and TBC when AM. In both cases, we define 
Nq = AM. When the system has an SU(2) symmetry relating to the spin freedom, the 
Fermi wave numbers for o =| and \ are same kp^ = kpi = hp, and we define the new 



numbers 28, BQ 



m 



= (N L ^ + N m + N Ll + N Rl )/2-N /2, 

n c = (N L ^-N m +N Ll -N Rl )/2, 

m s = (N L ^ + N m -N Ll -N Rl )/2, (27) 

= {N L ^-N m -N Ll +N Rl )/2. 



n 



2m c gives the total change of the electron number. The boundary condition of the phase 
field is given by 

(j) c , s {x + L) = (/) c , s (x) - V27rm C)S , d CjS (x+ L) = 9 C}S (x) - \/27m C)S . (28) 

Numbers in eq. fl2"7|) relate to the energy spectrum as 

(n-Lc, nuc, n Ls , n Rs ) - E 
2ttv c ( 1 2 i( c 2 \ 2ttv s / 1 2 K s 2 



2nv c , _, . 2tiv s _ 
— ^—{ricL + n cR )-\ ^—{n S L + n sR ), 



with the wave number 

( i 2tt \ 2tt 2tt. _ _ . . . 

q = - \ 2k F + -j-mc j n c - —m s n s - -j-( n cL - n cR + n sL - n sR ). (30) 

Relating to the SU(2) symmetry of the system, we have K s = 1. Integer numbers Nlv 
and N Rcr give the selection rule for the excitation: When 2m c is even (odd) integer, n c , 
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m s , andn s are integer (half odd integer). And when m c + m s is even (odd) integer, n c + n s 
is even (odd). 

The level spectroscopy with the selection rule has been applied for the spin-gap prob- 
lem of the t — J model [[29|, 0, and for the extended Hubbard model (half and quarter- 
filling) 0. 

For other applications, e.g., magnetic plateau, spin-Peierls transition, please see refer- 
ences in [K21 



6 Acknowledgement (en francais) 

Premierement, nous remercions Dr. K. Okamoto (application aux plateaux magnetiques) 
et M. Nakamura (application aux systemes electroniques) . 

Cette theorie etait nee, quand un des auteurs, K. N., etudiait en France, il y a sept ans. 
Le nom "level spectroscopy" etait nomme par professeur H. J. Schulz. Nous remercions 
profondement H. J. Schulz, parce qu'il a encourage K. N. a developper cette idee, et qu'il 
a apprecie cette theorie, apres 1'accomplissement. L'embryon de cette idee remonte aux 



recherches [031 IIT 



Indirectement, les discussions sur la correction logarithmique avec professeur I. Affleck 
par lettres (e-mail) ont aide le progres de notre recherche. Professeur M. Takahashi nous 
a invite a commencer l'etude sur la theorie des champs conforme. Les crititiques et 
les commentaires par professeur H. Shiba, par example, "C'est tres difficile a rechercher 
numeriquement la supraconduction, parce il y a la probleme de renormalization.", aussi 
sont de valeur. 
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